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Abstract 

In the conventional quantum mechanics (i.e., hermitian QM) the adia- 
batic theorem for systems subjected to time periodic fields holds only for 
bound systems and not for open ones (where ionization and dissociation 
take place) [D. W. Hone, R. Ketzmerik, and W. Kohn, Phys. Rev. A 
56, 4045 (1997)]. Here with the help of the (t,f) formalism combined 
with the complex scaling method we derive an adiabatic theorem for open 
systems and provide an analytical criteria for the validity of the adiabatic 
limit. The use of the complex scaling transformation plays a key role in 
our derivation. As a numerical example we apply the adiabatic theorem 
we derived to a ID model Hamiltonian of Xe atom which interacts with 
strong, monochromatic sine-square laser pulses. We show that the gener- 
ation of odd-order harmonics and the absence of hyper-Raman lines, even 
when the pulses are extremely short, can be explained with the help of 
the adiabatic theorem we derived. 

03.65.-w, 42.50.Hz, 42.65.-Ky, 32.80.Rm 



1 Motivation 

When matter is exposed to intense laser fields, high harmonics (HHs) of the 
incident radiation may be produced. Usually, only odd harmonics are obtained 
even when the laser pulses are short (for theoretical and experimental work 
which demonstrates this see £Q and |2] , respectively) . Since the duration of the 
pulse in time is inversely proportional to its width in energy space, one may find 
this result surprising, as one may expect to obtain also a large distribution of 
frequencies in the scattered field. Why are only odd harmonics obtained even 
when the laser pulses are short? 

For CW lasers (and symmetric field-free potential) using the non-hermitian 
Floquet theory it was proved that only odd harmonics are obtained when the 
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dynamics is controlled by a single resonance Floquet quasienergy (QE) state 
When laser pulses are used it was argued that this proof holds since usu- 
ally the populated resonance states are associated with very different lifetimes 
and the dynamics is controlled by the resonance state which has the longest 
lifetime. However, this argument may hold only when the duration of the laser 
pulses is large enough to enable the decay of the short lived resonances. Indeed 
numerical simulations showed that the harmonic generation spectra (HGS) as 
obtained from a single non-hermitian (complex scaled) resonance Floquet state 
is in a remarkable agreement with the results obtained from conventional (i.e., 
hermitian) time dependent simulation [5]. 

The question that is addressed in this work is weather an analytical criteria 
for the shape and duration of the laser pulse for which the system is controlled 
by a single resonance Floquet state can be given. It is obvious that the question 
regarding the possibility of population of a single resonance state is connected 
with the question regarding the degree of adiabaticity of the process. The 
question is therefore under which conditions can a short laser pulse be defined 
as adiabatic one. The answer to this question is important not only to harmonic 
generation (HG) studies but also for other, more general studies where lasers 
are used to control the dynamics, for example adiabatic STIRAP procedures 

In order to answer this question we use the (t,f) formalism [7] together 
with the non-hermitian quantum mechanics (NHQM) formalism. The use of 
NHQM formalism to describe the dynamics of atoms/molecules subjected to 
CW laser fields is essential, since only then can the dynamics be described in 
terms physical, square integrable, resonance Floquet states. Otherwise, the 
description of the dynamics in terms of hermitian Floquet states results in very 
little physical insight on the problem, as well as numerical problems, not to 
mention that it is limited to the description of bound systems only. In the 
hermitian case the spectrum is continuous and becomes discrete only due to 
the use of finite box quantization; moreover, a single Floquet state in hermitian 
QM can not describe neither the resonance phenomena nor the field ionization 
phenomena 8 . 

Our strategy is as follows. In section 2 we give a brief review of the for- 
malism used in our derivation [namely the (t,f) formalism for hermitian and 
non-hermitian Hamiltonians] . In section 3 we introduce the new derivation of 
the adiabatic theorem for open quantum systems in strong laser pulses. In sec- 
tion 4 we apply the adiabatic theorem as derived in section 3 to a test-case model 
hamiltonian which describes a ID Xe atom subjected to a sin-square pulse of 
monochromatic laser. In section 4 we conclude. 

2 A brief review of the (t,V) formalism 

The (t,f) formalism enables one to obtain analytical solutions for any time- 
dependent Schrodinger equation (TDSE) with time-dependent Hamiltonians. 
The formalism rests on lifting of the TDSE to an extended Hilbert space, prop- 
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agation of the wavefunction there and finally projecting back to the physical 
Hilbert space. 

The solution of a general TDSE is given as usual by: 



JI(r,t)tf(r,t) =ih— *(r,t) 



(1) 



with the initial condition ^(r, t = 0) and the vector operator r describes the 
internal degrees of freedom. It has been shown by Peskin and Moiseyev that by 
regarding time as an extra coordinate t' , one can obtain another Schrodinger 
equation with a time independent Hamiltonian in the extended Hilbert space 
(r, t'), whose solution *(r, t' , t) has an analytical time dependence (given by the 
analytical time evolution operator associated with time-independent Hamiltoni- 
ans) 0. Our desired wavfunction <3>(r, t) can be deduced from this wavfunction 
by a simple operation (that will be shown later). The advantage is that effi- 
cient propagation schemes designed for solving the TDSE with time- independent 
Hamiltonians could then be used also for time dependent Hamiltonians, thus 
releasing one from the difficulties associated with time ordering. 

Let us define the following Floquet-type operator 



H F (r,t') = (H(r,t)-ih-)\ t , =t = H(r,t')-ih— ,0<t'<t f (2) 



where t' should be regarded as a coordinate and tf is the final time of propaga- 
tion. When the Hamiltonian H (r, t) is time periodic with period T, the operator 
is the Floquet operator and tf =T. Provided that correct boundary conditions 
are chosen for the t' coordinate, this operator is hermitian since it is the sum of 
two hermitian operators. This operator has eigenstates and eigenvalues which 
are given by the eigenvalue equation: 



and the set of eigenstates is complete in the extended Hilbert space (r,i') with 
respect to the inner product 



Say the following TDSE with time-independent Hamiltonian need to be 
solved with the initial state ^(r, t' , t = 0): 



H F (r,t')(P a (r,t') = e a (f> a (r,t') 



(3) 




(4) 



H F (r,t')y(r,t',t)=ih—V(r,t',t) 
Using the definition of Hp(r, t') this equation reads 



(5) 




(G) 
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[H(r, 0* (r, t', t)]\ t , =t = H(r, t)[*(r, t)|f=t] = ^ ^ *)!*'=*] (7) 
where the following property 

|[W( r) t',f)| t , =t ] = [(| + |-)W(r,f',t)]| t , = , (8) 

(which holds true for any function oft' and t. The relation between the solution 
of Eq0and the solution of the original TDSE (EqQJ is given by 

tf(r,t)=¥(r,t',t)|t, =t (9) 

provided that these two differential equations have the same initial condition. 
Hence, 

*(r,t / ,t)| t , =t=0 = *(r,t = Q) (10) 

It is seen from Ea llOl that apparently the initial condition isn't unique. While 
this holds true in case that one is interested only in the physical wavefunction 
\I/(r,i), it should be noted that if one wishes to calculate physical quantities in 
the extended Hilbert space using the function *(r, t', t) and then go back to the 
original Hilbert space, and get the correct results, the initial condition in the 
extended Hilbert space should behaves as a delta- function in t' , The correct 
initial condition will therefore be *(r, if, t = 0) = *(r, t = 0)S(t'). 

The main advantage of the (t,f) formalism from a numerical point of view 
is that it enables the use of an analytical expression for the time evolution 
operator, without the necessity of time ordering, even when the Hamiltonian is 
strongly time-dependent. Any TDSE with time-dependent Hamiltonian could 
be replaced by a different TDSE, with time-independent Hamiltonian, for which 
a greater number of accurate integration schemes exist and the solution is given 
formally by ^'(r, t' , t) = e~^ HF ^ r,t )'\&(r, t', t = 0)). The price one pays however 
is that the new TDSE need to be integrated over one more dimension. 

The main advantage of the (t,t") formalism from a conceptional point of 
view is that it enables to describe any time-dependent dynamics in terms of 
stationary eigenstates and eigenvalues. 

So far the derivation has been carried out within the framework of the con- 
ventional (i.e., hermitian) quantum mechanics. Using the complex-scaling (CS) 
transformation |9l I10| 

H F (r,t') — ► H e F (r,t') = H(re ie ,t) - ih^ (11) 

the quasienergy spectrum of the Floquet Hamiltonian becomes complex and 
square-integrable resonance states, that were embedded in the continuum in 
the unsealed problem, are uncovered. For sake of simplicity we drop the index 9 
in all ^-dependent expressions (operators, eigenvalues, eigenvectors etc.) in the 
proceeding text. 
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The Floquet operator Hp can be represented with the orthogonal Fourier 
basis set -j=e lu)nt ,n = 0, ±1, ±2, ... as a square matrix 

i r T 

[ gf(r) W,n = ™ / dt'e-^ n 4 ff F (r, £')e" w (12) 
1 Jo 

where T = 2tt/uj. The left and right eigenvectors (bi-orthogonal set of eigen- 
vectors, of this Floquet matrix are the Fourier components of the Floquet 
states as defined in EqOU that is 

[H F {v)f^ L a {v) = Erffr) (13) 



PfWICW = E a l?K(r) (14) 

where 

W L J H {v)] n = <p L J* (r). (15) 

Since in our case [H F (r)Y = [Hp(r)] then y? a (r) — ~^a( r )- There are two 
sets of eigenfunctions of the Floquet operator Hp(r, t') 

^(r,t') = E<n( r ) eW ( 16 ) 

n 

and 

<^(r, = £ < n (r)e-^' = ]T < n (r) e --"*' (17) 

n n 

As pointed out in J2| the c-product, which is associated with the non usual 
inner-product in linear algebra (see for example Wilkinson's text book |ll|h 
reads 

(^(r,tO|^(r,t')) r ,t'=E«nl^,n) r = ^ (18) 

n 

where (^(r, OI^MOkt' = (r, Ol05( r >O)r,f • However, as proposed 
recently by Moiseyev and Lein (J3j), the inner- product should be modified 
even further when time-dependent functions are used as basis set due to the 
time-asymmetry problem in NHQM. If we define 0^ (r, t) and E a = E ra — 
to be the Floquet eigenfunctions and eigenvalues respectively, then the time- 
dependent basis functions could be 

*«(r,i^) = e-^V*(r,0 (19) 
(which are solutions of the TDSE) and also the left functions 

*L(r,t',t) = e + l E « t ri(r,t'). (20) 
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Following the modified definition of the inner product ("finite-range" product, 
"F-product" E3|) 

= e-* B - t e + * B :'*(^(r J Ol05(r,* / ))r,f = e^S^ (21) 

it implies that the a th quasienergy state decays exponentially in time. For more 
detailed discussion see [T3 ] .[T3 | . 

3 The adiabatic theorem for time-dependent sys- 
tems 

The adiabatic theorem for time-dependent bound systems was derived in 1997 
by Kohn et al. ^5] and in 1999 by Holthaus et al. JH] who used the (t,f) for- 
malism to describe the evolution of a system subjected to chirped laser pulses. 
As discussed by Kohn and co-workers, the adiabatic approach is not applica- 
ble for open systems since the quasienergy level spacing reduce to zero as the 
number of basis functions used in the numerical calculation is increased. To 
avoid this difficulty Baer et al. ^7] applied the adiabatic theorem to time- 
dependent open systems in the high frequency regime where the system was 
stabilized and the resonances (which were embedded in the continuum part of 
the Floquet spectra) became so narrow that they could be practically treated 
as bound states. The purpose of our work is to derive the adiabatic theorem for 
general time-dependent open systems where there are no bound states and the 
resonances aren't necessarily narrow. We are using the non-hermitian Floquet 
formalism (through the CS formalism) which allow us to describe the dynamics 
in term of non-hermitian resonance states (see for example [9"|-|21|. and also the 
work of Day et al. |22 who used the NH Floquet multistate method to study the 
applicability of the single Floquet resonance approximation in the description 
of the dynamics of H atom subjected to intense laser fields of various strengths). 

Below we derive the adiabatic theorem for time-dependent open systems 
using the extended (t,f) formalism. By the term "extended" (t,f) formalism 
we mean that in the same manner presented, one may add any number of time 
"coordinates" to the Schrodinger equation as one wishes, if by this a better 
understanding or easier solution of the problem is achieved. Here we found that 
by addition of 2 time " coordiantes" to the TDSE, we simplified the derivation 
of the adiabaticity criteria for photo-induced dynamical systems. In this sense, 
we are using a (t,t',t") formalism. 

We would like to study the dynamics of a single active electron in an atom or 
molecule, subjected to a pulse of strong monochromatic linearly polarized laser 
radiation. In the dipole approximation the TDSE which describes this process 
is: 

H(v,t)V R (r,t) =ih^ R (v,t) (22) 
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where, 



H(r, t) = H (r) - er • f(t)cos(ut) (23) 

and, 

f (t) = e ek/(t). (24) 

Here f(t) is the function which describes the envelope of the laser pulse and 
f(t) is a vector as defined in Eq|22| e is the laser's amplitude, ek is a unit 
vector in the direction of the electric component of the laser field, li is the 
laser's frequency, with T = 2tt/oj the optical period. H (r) is the field-free 
Hamiltonian and the vector operator r describes the internal degrees of freedom 
(the coordinates are complex-scaled throughout) . 

In the same spirit of subsection (2.1), we define the following operator 

H F (r, t', f) = H(r, t', t") - ih^ - ih-^ (25) 

where 

H(r, t', t") = H (r) - er ■ i{t")cos{ut') (26) 

and t', t" should be regarded as additional coordinates. Upon complex-scaling 
|2*3~] Hp(r,t' ,t") becomes non-hermitian. Therefore the inner c-product should 
be used as mentioned before. The quasi-energy solutions of this complex-scaled 
Floquet-type operator are: 

H F (r,t',t")^(r,t',t") = X k ^(r,t',t") (27) 

H F *(r,t'X)^(r,t'X) = \ k ^(r,t'X) (28) 

where the symbol Hp doesn't stand for an operator but for the transpose of 
the matrix representing the operator Hp. The eigenfunctions form a complete 
set in the extended Hilbert space r,t',t". 

Say we want to solve the following TDSE with the initial state ^'(r, t' ,t" ,t = 

0): 

H F (r, t', i")$ fl (r, t', t", t) = ihj* R (r, t', t" ,t) (29) 
The solution to this equation is 

$ fl (r, t', t", t) = e -kHH*,t',t")tyR^ ^ ^ t = 0)=Y, c fc e-* Afe Vf(r, t', t") 

k 

(30) 

and a function fy L (r, t', t", t), which is not a solution of any Schrodinger equa- 
tion, is defined as 



7 



^ L (r,t',t",t) = J2^e + ^ t ^(r,t'X) (31) 

k 

where by taking the cut t' = t" = t = on Eq[3U| it is easily seen that the 
expansion coefficients are c k = (ip% (r, t' , t")| 4 / =4 » =0 |* R (r, t' , t" , t)\ t > =t "=t=o)r- 
[e~i Xkt ip£-(r,t',t")\t>=t»=t is a solution of the original TDSE (apply the cut 
t'=t"=t on Eg 1271 and compare the result to the result obtained when the func- 
tion e~Ti Xkt ?pJ*(r, t', t")\ t '=t"=t is substituted in the original TDSE); therefore, 
any linear combination of these solutions is also a solution] . 

Let us now return to the main purpose of this article, the derivation of the 
adiabatic theorem for non-hermitian open systems. We would like to treat t" as 
an adiabatic coordinate (this is the " coordinate" associated with the pulse enve- 
lope) in the same way that the electronic motion is separated from the nuclear 
one in the treatment of molecules within the Born-Oppenheimer approximation. 
First we define the following operator 

H ad (r,t',t")=H(r,t'X)-iK-^ (32) 

where t" should be regarded as a parameter now. This means that this Hamilto- 
nian is a Floquet Hamiltonian describing the interaction of the atom with CW 
laser of strength e\ where following Eal24l 

e x = |f(t)| = eo/(t) (33) 

The eigenstates of this operator form a complete basis (in the r — t' space), 
for every value of the parameter t" : 

H ad (r, t>, t")r a d ' R (r, t>, t") = e a a d (t"W a d ' R (r, f, t"). (34) 

Notice that due to the complex scaling, Afc and s% d (t") get complex values. 
We can expand each eigenstate of the complete problem (Eg 1270 in this basis: 

<(r, ;',;") =J2rJ' R (r,t',t") Xa 'Mn (35) 

a ' 

Substituting Eq|23into EqE| multiplying the obtained equation from the left 
hand side by ifi^,' (r, t',t") and integrating over r and t' one gets, in matrix 
notation, the equality: 

l-i^L+ (E ad {t")+Yit")Mk{t") = Afc* fc (t") (36) 

where 

{E ad (t")} a ,a' =e a J(t")6 a , a ,, 
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[**(*")]» =X«,*(*")- (37) 

Notice that in case that the matrix on the left hand side of Eq|2S]is diagonal, 
a homogeneous systems of uncoupled equations is obtained. In such a case one 
should solve each equation separately. Therefore, the sum in Eg 1351 reduces to 
a single product. This is exactly the adiabatic approximation as appears in the 
Born-Oppenheimer context. The next step in our derivation is to represent the 
matrix (E_ ad (t") +V_(t")) by its spectral decomposition 

\E ad (t") + Yit")}R R {t") = Il R {t")Mt") (38) 



\E ad (t") +YSt")} % L {t") = nr(t")K(t") (39) 

The matrix of eigenvalues W(t") is diagonal and the right and left eigenvectors 
are normalized with respect to each other in order to maintain the correct inner 
product: 

[R L (t")] t R R (t")=l (40) 

In the case that the matrix (E_ ad (t") + V_(t")) is not strictly diagonal we can 
use first-order perturbation theory to get the first-order deviation from diagonal. 
If we treat the matrix V_(t") as perturbation, we get 

The matrix will be diagonal to a good approximation if 

A a (t")= IL£V')W,a|«l (42) 
which produces the following adiabaticity criteria, 

A (f „, v ((cWr)|-^IC?M^0))^ „ . ( ^ 

M } = & 1 1 <K (43) 

Using the specific form of the Hamiltonian of the problem given in Eg 1231 
and the Hellman-Feynman theorem one gets the adiabatic condition for time- 
dependent open systems: 

A Q (t")«l (44) 

where 



A «»\ I \%, I #(*"), \- («' L (r, t', t")\r ■ e k cosM')|C^(r, t t"))) rA , 

~& r Km-sfit")? 1 

(45) 
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The index a is a super-index; since fEa l34f> it is easily seen that not only is 
ip%?' L ^ R (T,t' ,t") a solution of the eigenvalue equation, with eigenvalue e^ d (t"), 
but also e lujmt ' ipa d ' L ^ R {v, t', t") is a solution, with the eigenvalue e^ d (t") + htum, 
for any integer m. Let us take all the states whose corresponding eigenvalues lie 
in the interval [0, tko] (the first Brillouin zone) and define them to have index 
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0; we call these states ' (r,t',i"), 0* : (r,t',t") and the corresponding 
eigenvalues E°- d (t") and get: 

Tp° d < R (r, t\ t") ee (/>f' R (r, t', t ")e lujmt ' (46) 

C d ' L (r, t', t") ee <f>f' L (r, t', t")e™ mt ' (47) 

e ad {t") ee ^ Qrf (t") + Sum (48) 

where 

< Ef(t") < Tilu (49) 

Thus, the index a actually counts both the position of the quasienergy within 
the first Brillouin zone (the index j) and the Brillouin zone itself (the index n). 
With respect to the generalized inner product, two states with one or more of 
the indices (j,n) different arc orthogonal. 

Going back to Eg 1451 now ,it is seen that the probability to couple an initial 
adiabatic state ip^Ltj j (r, t', t") to any other adiabatic state ip^_>; m ,\ (r, t', t") 
is given by, 



A (m (t")=F(t")-^a^ ) (t") (50) 



where 



F(t") = \e\ns \^-\, (51) 

«$)(0 =E (52) 

and the functions Cy, ^(t"), that will be termed as "adiabatic cross terms" 
from now on, are given by 

m ((<^^r,t^^0|r^e k c OS (a;Oe^>° d ^ fl (r,^,t"))) ^ , t , 

j -Ef{t") +Humi) 2 1 J 

Since the energies E°; d (t") are complex (the Hamiltonian is non-hermitian) 
then for j ^ j' it is most unlikely that E°", d = E^ d . It is clear that the denomi- 
nator hardly ever vanishes even when m = 0. This holds true even when j = j' 
but m ^ 0. 
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Therefore, the criteria for a pulse to be considered adiabatic is that the con- 
dition 

A (j . 0) (O«l (54) 

be fulfilled. 

Who is the adiabatic state V'af^j o)( r > whose couplings to all other 

states m n(r, t' , t") should remain small in the adiabatic limit? Assuming 

that the system is in a stationary state of the field- free problem <pj (r) before 
the action of the field (relaying on the superposition principle of the solutions 
of the TDSE generality is not lost by this assumption) and provided that the 
field is switched adiabatically, this Floquet resonanace state tp^_/j (r, t\ i") is 
the state which is "born" from the stationary state fj{r) as the field is switched 
on. If the process is not done adiabatically many Floquet resonanace states will 
be populated, resulting in considerable couplings of iplfLu )( r i t\ t") to them 
and collapse of the adiabatic condition in EqEH The only adiabatic check 
which is physically meaningful is one in which a denotes a resonance state 
(which is associated with square-integrable function), a' however stands for 
both resonances and rotated continuum states. 

One should notice that the derivation of the adiabatic theorem presented 
above holds for many electron systems. In order to avoid complicated notation 
the symbol r can stand for many electrons. It also holds for polychromatic radia- 
tion, where the CW field is a collection of monochromatic fields with frequencies 
oji and phases tfi. 

Notice that for a given problem (given spectral profile of the CW field and 
given field- free potential), the sum over absolute value of the adiabatic cross 
terms should be calculated as function of the effective CW-field intensity (which 
is symbolized here through t") only once. The adiabatic cross terms should 
then be converted to be functions of time, through the explicit time dependence 
of the pulse envelope and then the sum of their absolute values should be mul- 
tiplied by the time derivative of the pulse envelope and by the maximal field 
intensity to obtain the final expression which indicates whether the adiabatic 
criteria is fulfilled or not. 

It is easily seen in Eg 1511 that for a given system, the shape and intensity of 
the laser pulse determines its adiabaticity since these parameters influence the 
shape-derivative term. A short pulse which is switched on or off abruptly and 
has a high maximal intensity, will most likely not be adiabatic. 

In the case that the adiabaticity criteria is fulfilled, the sum in Ea 1351 could 
be reduced to a single product: 

^ 3 (r,t'X)^^f R (r,t'X) Xj Mt") (55) 

and the adiabatic states are assigned with two good quantum numbers k and j. 
The solution to the eigenvalue equation Eq|3fj|is now {\Vit")] a . a ' ~ 0) 

Xjtk (t") = e+i x >° t "e-xr"< 1 -<V. (56) 
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By using Ea l3QI55l56l we get that the adiabatic solution, in the (t,t',t") formal- 
ism, is given by, 

y R (r,t',t",t) = Ecje-W^')]^ ' R (r,t',t")e-%f " drE ? d( - T) (57) 

k 

Now, applying the cut t" — t in order to obtain the physical solution of Eal22l 
one gets (eliminating the phase factor ^ fc c^): 

**(r,t) = * R (r,t',t",t)\ t „ =t , =f = e-*r d - E ° d ^cbf R (r,t',t")\ t „ =t ,= t (58) 

This is the adiabatic solution of the TDSE associated with initial state (pj(r). 

Let us summarize and clarify the procedure that need to be made in order 
to determine if a pulse is adiabatic or not; The determination is carried out 
through the calculation of the expression Atj ^(t) (it is t now, not t" !): 

(1) Perform non-hermitian adiabatic Floquet simulations fEa l34|) with CW 
field, for a range of intensities £\ which covers all intensities between zero and 
the maximal intensity Eq that the studied laser pulse reaches. The adiabatic 
Floquet Hamiltonian is therefore 

Q 

H ad {v, t', ex) = H (r) - er • Eicc>s(u)t ) ~ ifi— (59) 
and the eigenvalue equation to be solved is 

H 0< ,(r,* , ,e 1 )^' Jl (r.* , ,ei) = Ef(E 1 )(bf R (r,t' 1 E 1 ) (60) 
Obtain the quasienergy spectrum Ej d (£i) and all the adiabatic cross terms 

c (j'°m)( £l ) as defined in Eg 1531601 as function of the intensities. This stage is 
done only once, for a given system. 

(2) For a given laser pulse f(t) with a maximal intensity Eq, evaluate the 
effective CW-field intensity as function of time £of{t) fEa 124(1 . Then, convert 
the adiabatic cross terms to be functions of time via the equality (EqESJl 

£i = £of{t)- 

using c^',°^(ei) as calculated in step 1 c^°^(<) = c[i'°m)(/ _1 ( £ i/ e o)) wnere 
/ is the transformation which fulfills / -1 [/(i)] = t. 

(3) For each given resonance state a = (j, 0) calculate 

%o)W^W'EEiS)Wi (6i) 

using c|^',°^(i) from step 2 where notice that here t" in Eas. (|5QI53|l is replaced 
by t. If for a given resonance state a — (j, 0) the corresponding expression 
A(j z o)(t) is smaller then unity at every instant, it is guaranteed that the system 
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initially at the bound state which corresponds to this resonance, will evolve 
adiabatically to that resonance. In this case the HGS spectra will show only 
odd harmonics and the ionization probability as function of time will have a 
simple form that will be shown. 

4 Illustrative numerical example 

We studied a single-electron ID Xe atom subjected to a single sin-square pulse 
of strong monochromatic laser field in two approaches. In the first one hermitian 
simulations were carried out whereas in the second non-hermitian Floquet sim- 
ulations based on the complex scaling method were carried out. The hermitian 
simulations were carried out by solving the following TDSE 

[-—g-2+V (x)-exE f(t)cos(cjt)Mx,t)=ih^(x,t) (62) 
with a sine-square envelope 

/(<)=sin O (63) 

and the field- free effective potential Vq(x) was an inverse Gaussian 

V (x) = -0.63exp(-0.i424iE 2 ) (64) 

which supports two bound states that mimic the two lowest electronic states 
of Xe, with energies Eq = — 0. 4451a.it., E\ — — 0.1400a.w. and a third weakly 
bound state with energy E2 = — 0.00014a.u. 

The wave function was taken initially at the ground state (g.s.) of the field- 
free Hamiltonian of the system: 



V(x,t = 0)=<p 1 (x) (65) 

and was calculated for times < t < NT (single sine-square pulse, N was the 
number of optical cycles that entered the pulse). Using this wavefunction the 
ionization probability at times t = ^ + n^,0 < n < 2N — 1 was calculated (at 
these times the potential felt by the electron was the field- free potential Vq(x)) 

3 

p ion (t) = i-Y / \(<Pi\nt))\ 2 (66) 

i=l 

where ipi were the bound states of the field-free problem (3 bound states over 
which we summed in this example). Also HGS was calculated, which follow- 
ing the classical-quantum correspondence principle (Larmor formula. |24J) equals 
the modulus-square of the Fourier-transformed time-dependent acceleration ex- 
pectation value. This is actually the intensity of the radiation emitted by the 
oscillating electron as presented in energy space. 
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a(n) = \—j o -^(t)\--^ + ee f(t)cos(cJt)m)) e~ mt dt\ 2 (67) 

The non-hermitian simulations were Floquet simulations which were carried 
out for different field intensities. The quasienergy spectrum of complex energies 
E°-, d (ei) and the adiabatic cross terms c ui°ln){ £ i) were calculated for each inten- 
sity. Then, these quantities were expressed as function of time through Eg 1331 
i.e. t = / _1 (ei/eo)- It was verified that in the cases where the adiabatic crite- 
ria was fulfilled, the HGS obtained from the hermitian propagation simulation 
contained only odd harmonics. A more quantitative measure of the existence of 
the adiabatic criteria was obtained through the comparison of ionization prob- 
abilities as obtained from the hermitian propagation simulation and the non 
hermitian simulation. 

In order to get non-hermitian Floquet Hamiltonian, the complex coordinate 
method was used. The Floquet Hamiltonian (EqEHJ 



fl od (a; 1 t , ,ei) = '^^q^ + V (xe a ) - exe i6 e f(t)cos{u>t) - ih— (68) 

was diagonalized. Provided that the scaling parameter 6 was sufficiently large, 
the resonance quasienergy states were 6 independent: 

Ef{t")=E ri (i")-jli(t") ; T j = - (69) 

Z Tj 

E r . being the position of the state, Tj being the width of the state and Tj its 
lifetime. Since the resonance states had finite lifetimes, and since the resonances 
are the states which are associated with the dynamics, these resonance lifetimes 
should have fingerprints in the hermitian propagation simulation. Indeed it 
was found to be so when the ionization probabilities, as computed in the two 
simulations, were compared; The ionization probability at each instant was given 
by the following expression, which was obtained using Eg 1581 and the F-product 
definition for the inner product |13|: 

Pion(t)=e-*f dTr ^ T) (70) 

where T\(r) was associated with the Floquet resonanace state that was "born" 
from the ground stationary state <fi{v) as the field was turned on. The res- 
onances that were "born" from the field- free Hamiltonian bound states were 
identified by plotting the quasienergy spectrum as function of the effective field 
intensity £\. The resonance complex quasienergy trajectories started from the 
field-free Hamiltonian bound states real energies and formed continuous trajec- 
tories in the complex energy plane as function of the intensity. 

How was this expression for the ionization probability obtained? Accord- 
ing to the F-product definition when the complex energy given in Eg 1691 is 
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substituted in Eq|5S| the "ket" (right) solution of the non-hermitian TDSE is 
obtained: 

V R (r,t) = e -^f dTT i {T) -*/ t,lrB 'i (T) 0° < '" B (r,t',t")|t«=t'=f (71) 

and it is easily seen that this function decays with time. The "bra" (left) 
solution is not a solution of a Schrodinger equation but is derived from the 
"ket" solution (the explanation of this point is beyond the scope of this work; 
for an explanation see |14p 

* L (r,<) = e -*f dTT ^ T) + *f tdTEr > (T Uf' L (r,t',t")\ t „=t>=t (72) 

and also this function decays with time. When the F-product of ^ L (r,t) and 
^> R (r, t) is calculated (this is an overlap integral without complex conjugation of 
the left state; in the usual dirac-product notation this reads (^ L *(t)\^ R (t))) the 
terms containing the real part of the energy cancel each other and the overlap 
integral (c-product) of the adiabatic Floquet states gives unity (remember the 
completeness property of Floquet states also in coordinate space alone). The 
expression in Eq[7Ulis obtained for the specific case that only the Floquet reso- 
nance state which is "born" from the field- free bound state is populated. When 
several Floquet states are populated the same type of calculation could be re- 
peated, where this time ty L (r,t) and ty R (r,i) are given by liner combination 
of terms as in EqEJwith the proper coefficients. Since the adiabatic Floquet 
states are orthogonal to each other, a generalization of the result of Ea l7UI is 
obtained: 

3 

JW*)=£l<¥>i|*(t = 0))| a e-*J dTTi ^ (73) 
i=i 

Notice that in the non-hermitian formalism the norm is not conserved but decays 
with time; it contains the knowledge about decay processes inherently. 

The two functions as given in Eg 16617(71 were compared and it was found 
that the resemblance between the functions increased as the adiabatic limit was 
increasingly reached by the laser pulse parameters. 

The numerical method used to solve the TDSE with hermitian Hamiltonian 
was the split operator Forest-Ruth algorithm with 7 points |25| : The grid size, 
time step and/or grid step were adjusted as required to achieve convergence. 
For the non-hermitian Floquet simulation, the (t,f) formalism was used, with 
the complex coordinate method. The number of basis functions, box length and 
scaling angle were adjusted as required to achieve convergence. 

In Fig^the HGS as obtained for pulse strength of Eq = 0.035a. it. (corre- 
sponding to intensity of 4.30 -10 13 PF/ cm 2 ), pulse-durations of N = 5, 10, 50 opti- 
cal cycles, laser frequencies of u> = 0.015a. it., 0.07a. it., O.lla.u. (corresponding to 
energies of Ticu = 0.408eV, 1.90eV, 2. 99eV and wavelengths of A = 3040nm, 651nm, 
414nm, respectively) and initial state <pi(x) are shown. The appearance of odd 
harmonics and the absence of even or non-integer harmonics, is in general the 
main feature which appears for long pulses, regardless of the frequency. The odd 
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harmonics arc obtained even for the shortest laser pulses although the existence 
of odd-symmetry selection rules is sensitive to the frequency: some deviation 
appears at frequency of lj = 0.07a. it. and the obtained spectra is more compli- 
cated. 

In Fig [21 the complex quasienergies of the 3 resonances which are "born" 
from the 3 bound states of the field-free Hamiltonian (will be given indices 1- 
3 from now on), as obtained from adiabatic Floquet simulations fEa l60|) for 
lj = 0.11a. it. are shown as function of the CW-field strength E\. Notice that the 
lifetimes of these resonances are not monotonic functions of the field intensity. 
It should be noted that these were not the only resonances that appeared in the 
quasienergy spectrum; there were also other resonances which didn't emerge 
from the field-free bound states. However, as will be shown in FigQJ in this 
case for not too large field intensities resonances 1-3 were the only important 
resonances and they alone determined the dynamical behavior of the system. 

In FigElthe ionization probability as obtained from the hermitian simulation 
fEa l66[l with pulse strength of Eq = 0.035a.u., pulse-durations of N = 5, 10, 50 
optical cycles, laser frequency of lj = O.lla.u. and initial state ^1(2;), is com- 
pared to the ionization probability as obtained from the expression given in 
Eg 1701 which is derived from the F-product formalism together with the res- 
onances quasienergies obtained from the non-hermitian simulation. It is seen 
that as the pulse becomes longer the results obtained from the two simulations 
become identical. 

Both results of HGS and ionization probabilities showed that the time- 
dependent wavefunction of the studied systems could be well approximated by 
the adiabatic expression given in EqEHI even for short pulses. The values of the 
terms which check the adiabatic criteria, as seen in FigEl gave the explanation 
why this was so. 

In Fig0] the expressions a^j ; a (3) 1 a (3) (EqEJ which describe the couplings 
between every two resonances from the set of 3 tracked resonances is shown 
as function of the CW-field strength ex, for the case ui = O.lla.u.. It can be 
seen that the coupling between resonances 1 and 2 is strong for field strength 
of £1 ~ 0.04a. it., 0.064a. it. and 0.072a.u.. This could be partially explained on 
the basis of the values of the quasienergies, as seen in Fig0 at least for two 
field strengths out of the three. It is seen that for field strength of 0.04a.w. and 
0.064a.u. the real parts of the quasienergies cross, resulting in small value of 
the denominator in the expression given in Ea 1531 for cS'^-j (at least for one m 
term). In the same way the strong couplings between resonances 1 and 3 at field 
strength of e\ = 0.064a.u., and between resonances 2 and 3 at field strength of 
ex = 0.066a.u. could be explained on the basis of the quasienergy values at these 
field-strengths. In particular it should be noticed that the couplings between 
resonances 1 and 3 at ex = 0.064a. it. are the strongest among all 3 resonances 
due to the close values of both real and imaginary parts of the quasienergies. 
However, it should be noted that crossings in the quaisenergy plot are not always 
indications of large couplings since also the overlap between the wavefunctions 
(the nominator of the expression for c^f^s) is important. 
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In FigOthe sum Ylj'jti a (j')( £l ) wmcn describes the couplings between all 
quasistates of the system a' — (j' , m) to resonance 1 a = (1,0) is shown as func- 
tion of the CW- field strength e±, for the case uj = O.lla.u.. In addition, also the 
partial sum a^j (ei) + aSj (ei) which describes the couplings between resonances 
2 and 3 to the resonance state a = (1, 0) is shown. It can be shown that up to a 
moderate intensity of e\ ~ 0.04a. it. the first resonance is mainly coupled only to 
the other 2 resonances and not to other, higher resonances or continuum states. 
The entire dynamics is governed almost solely by the 3 resonances which are 
"born" from the 3 bound states of the field- free Hamiltonian. 

The sums described above, which are functions of a CW-field strength e\, 
are converted to be explicit functions of time for the specific sine-square pulse 
with maximal intensity Sq = 0.035a.u. used in the simulation (FigEI upper 
part). For this purpose for each time t < t < NT the effective CW-field 
strength 0.035/(i) is calculated and the values of the 2 functions shown in Fig0 
which fit this effective CW-field strength are taken. Hence the 2 sum functions 
are converted to be explicit functions of time. In the middle part of FigEI the 
full term Anm(t), which represents the degree of adiabaticity in the process of 
shining a ID Xe atom initially at the ground state with a sine-square laser pulse 
supporting N — 5 optical cycles of monochromatic radiation with frequency 
uj = 0.11a. it. and strength Eq = 0.035a. it., is shown as function of time. It is 
seen that the term A(im(t) is bounded by the value of ~ 2 • 10~ 3 for all times, 
whether it is calculated by coupling of the first resonance to other 2 resonances 
only or to all other states. The structure of the HHS for N = 5 seen in FigQ 
implies that this value is indeed small and the process is adiabatic. In the lower 
part of FigOthe same quantity is shown, but for to — 0. 07a.u.. Here it is seen the 
terms A^ 10 ^(t) are bounded by a much larger value of ~ 10 _1 for all times and 
the more complex structure of the HHS for N — 5 seen in Fig^inrplies that this 
value is not small enough to indicate the appearance of an adiabatic process. It 
should be noted that for a given system with field-strength £q, frequency uj and 
sine-square pulse envelope, we have F(t) — \e\(^f^-)\sin(j^)\. Therefore, as the 
number of optical cycles the pulse supports increases the general shape of the 
terms ^4(i,o) (*) is kept the same but is attenuated. Therefore, for uj = 0.07a.u. as 
N gets bigger, the pulse's envelope varies more slowly and the process becomes 
more and more adiabatic, as seen in the HHS spectra for N = 50 for example. 

5 Conclusions 

With the help of the (t,f) formalism we derive an adiabatic theorem for open 
systems. The use of the complex scaling transformation plays a key role in 
our derivation. For example, the spectrum of the Floquet Hamiltonian of 
an open system is changed dramatically. Rather than a continuous spectrum 
that is responsible for the absence of an adiabatic limit for TV (number of basis 
functions)— > oo in the conventional QM, the resonances are associated with a 
point spectrum and are separated from the continuum which is rotated into the 
lower half of the complex energy plane. 
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An interesting important numerical result of our derivation is that the cal- 
culation of the effect of the pulses's shape on the dynamics does not require 
heavy computations. The entire effect of the laser pulse is embedded in a mul- 
tiplication factor of df(t)/dt where £of(t) is the variation of the maximum field 
amplitude as function of time. 

As a numerical example we applied the adiabatic theorem we derived to 
a model Hamiltonian of Xe atom (with symmetric field-free potential) which 
interacts with strong, monochromatic laser pulses. We have shown that the 
generation of odd-order harmonics and the absence of even-order harmonics, 
even when the pulses are extremely short, can be explained with the help of the 
adiabatic theorem we derived. 

The use of a single-electron ID model to describe a realistic atom is justified 
since it has been shown before in many cases that all the main strong field 
effects are reproduced. Therefore the conclusions obtained with this model are 
also valid for a realistic atom. 
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Figure 1: (color online) HGS obtained from hermitian simulations describing 
a ID Xe atom subjected to various sin-square pulses with strength of £o = 
0.035a.u., laser frequencies of u> = 0. 015a.M. (upper part), u> = 0.07a. it. (middle 
part) and oj = O.lla.u. (lower part) and pulse-durations of TV = 5 (solid), 10 
(dot-dashed) and 50 (solid-dotted) optical cycles. The initial state was taken to 
be the ground state in all simulations. For the longest pulse duration only odd 
harmonics appear in the spectrum, and hyper Raman lines are absent. This is 
also the basic feature in the case of short pulses, although some deviation from 
this structure appear at frequency of w = 0. 07a.u.. 
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Figure 2: (color online) Positions (upper part) and lifetimes (lower part) of the 
complex quasienergies of resonances 1-3 as function of the CW-field strength 
£\ as obtained from adiabatic non-hermitian Floquet simulations fEa l60l6§|l for 
u) = O.lla.u.. 
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Figure 3: (color online) Ionization probabilities obtained from hermitian sim- 
ulations (lines) (EqJHBJ) for pulse-durations of N = 5 (upper part), N = 10 
(middle part), N = 50 (lower part) optical cycles, and from non-hermitian sim- 
ulations (dots) (Eq[7U|). As the pulse becomes longer the results obtained from 
the two simulations become identical since the adiabaticity of the process is 
increased. 
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Figure 4: (color online) The expressions 6o| 2 j(ei) (solid), aLj(£i) (dashed) , 

2.5a|gj(ei) (dotted) fEa!521) as function of the CW-field strength s\, for uj — 
O.lla.u.. 
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Figure 5: (color online) The expressions J2j'^i a y/) ( £ i) (dashed) and a^j (e±) + 
a (3)( £l ) (solid), which describe respectively the couplings between all quasis- 
tates of the system a' — (j',m) or only resonances 2 and 3 to resonance 1 
a = (1,0) are shown as function of the CW-field strength e\ 1 for the case 
uj = O.lla.u.. 
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Figure 6: (color online) upper part: The time-dependence of a sine-square 
pulse Eosin 2 (j£j) (solid) and a scaled time derivative 90(^)sm(^) (dashed) 
for N — 5, lu = 0.11a. it. and Eq = 0.035a. it.. Middle part: the full term 
■4(1,0) (*) (solid), and the partial term F(t)(a^(t) + a9^(t) (dashed) are shown 
as function of time. Both terms are bounded by the small value of ~ 2 • 10~ 3 
for all times. Lower part: the same as in the middle part, but for lu = 0.07a.u.. 
Both terms are bounded by a larger value of ~ 10 _1 . It is therefore deduced 
that for lu = O.lla.M. the process is adiabatic but for lu = 0.07a.u. it is not and 
this is indeed verified in the HGS given in FigQ]for N = 5. 
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